clc;
clear;

% 定义目标函数  
myfun = @(x) -x(1)* x(2)*x(3);  
% 定义约束条件  
A = [1,2,2;-1,-2,-2];  % 系数矩阵  
b = [72;0];    % 约束条件右端  
% 定义起始点  
x = [10, 10, 10]; 
% 求解优化问题  
[x, y] = fmincon(myfun, x, A, b, [], [], [], []);  
% 打印结果  
disp('X的值');disp(x);  
disp('函数值');disp(y);

% 使用CVX进行求解
x1 = 10;
x2 = 10;
x3 = 10; % 起始点
cvx_begin
    variables x1 x2 x3 ;
    maximize (log(x1) + log(x2) + log(x3));
    subject to
        x1 + 2.*x2 + 2.*x3 <= 72;
        x1+2*x2+2*x3>=0;
cvx_end
fval=-exp(cvx_optval)